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Abstract. We calculate running coupling corrections for the lowest-order gluon 
production cross section in high energy hadronic and nuclear scattering using the BLM 
scale-setting prescription. At leading order there are three powers of fixed coupling; 
in our final answer, these three couplings are replaced by seven factors of running 
coupling: five in the numerator and two in the denominator, forming a 'septumvirate' 
of running couplings, analogous to the 'triumvirate' of running couplings found earlier 
for the small- rr. BFKL/BK/JIMWLK evolution equations. It is interesting to note that 
the two running couplings in the denominator of the 'septumvirate' run with complex- 
valued momentum scales, which are complex conjugates of each other, such that the 
production cross section is indeed real. We use our lowest-order result to conjecture 
how running coupling corrections may enter the full fixed-coupling fcr-factorization 
formula for gluon production which includes non-linear small- a; evolution. 



1. Introduction 

This proceedings contribution is based on pQ. 

While an exact analytic formula for gluon production in AA collisions in the Color 
Glass Condensate (CGC) framework is still not known, we do know that in pp and 
pA collisions gluon production at the level of classical gluon fields and leading-ln 1 fx 
nonlinear quantum evolution is given by the kr- factorization formula |2, 3J: 

Here Y is the total rapidity interval of the collision, Cp = {N 2 — 1)/2N C , boldface 
variables denote two-component transverse plane vectors k = (A; 1 , A; 2 ), and <f> p , <$a are 
the unintegrated gluon distributions in the proton and the nucleus, respectively, which 
are defined by [2] 

Mk,y) = — I d 2 bd 2 re-* kr V 2 N G (r,b,y) (2) 
a s (Air) 6 J 

and 

<f> p (k,y) = — ^ / d 2 bd 2 re~ lh - r V 2 n G (r,b jy ). (3) 
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In Eq. (J2J) the quantity N G (r, b, y) denotes the forward scattering amplitude for a gluon 
dipole of transverse size r with its center located at the impact parameter b scattering 
on a target nucleus with total rapidity interval y. iV G (r, b, y) can, in general, be found 
from the JIMWLK evolution equation. In the large- N c limit it is related to the quark 
dipole forward scattering amplitude on the same nucleus N(r, b, y) by 

N G (r, b,y) = 2 N(r, b, y) - N(r, b, y)\ (4) 

where N(r, b, y) can be found from the BK evolution equation. The quantity na(r, b, y) 
from Eq. ((3]) is also a gluon dipole amplitude, but taken in a dilute regime, where it is 
found by solving the linear Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equation. 

Eq. ([I]) for the gluon production was derived in the fixed coupling approximation. 
However, the dipole amplitudes N(r,b,y), Na(r,b,y) and nc(r,b,y) are now known 
for the running coupling case due to the completion of the running coupling calculations 
for the BFKL/BK/ JIMWLK evolution equations (4j|5]. Using the running-coupling BK 
(rcBK) equation, the calculational framework presented above, when applied to heavy 
ion collisions by replacing p — > 0^ in Eq. (JTJ, and implemented with a careful inclusion 
of the nuclear geometry fluctuations, led to the prediction made in [6] of the charged 
particle multiplicity as a function of the collision centrality for the LHC This prediction 
was confirmed by the ALICE data in [7]. However, at the time the prediction [6] was 
made, the scales of the couplings explicitly shown in Eqs. (JTJ, (|2]), and (J3j) were not 
known and had to be modeled. Our goal here is to fix the scales of those couplings for 
the lowest-order gluon production. 



2. Running coupling corrections: strategy 

To include running coupling corrections we follow the BLM scale-setting procedure [8] , 
which is known to be correct at least at the leading order in 1/N%. One first needs 
to resum the contribution of all quark bubble corrections giving powers of a^Nf, 
with Nf the number of quark flavors and the physical coupling at some arbitrary 
renormalization scale \i. We then complete Nf to the full beta-function by replacing 

N f ->-6irp2 (5) 
in the obtained expression. Here 

ft = (6) 

is the one-loop QCD beta-function. After this, the powers of @2 should combine into 
physical running couplings 

asiQ2) = rz i q 2 (7) 

1 + a M 02 In ^2- 

at various momentum scales Q which would follow from this calculation. We use the 
MS renormalization scheme. 

As was originally argued in [9], including running coupling corrections into the 
diagrams of Fig. [TJ below assuming that only a gluon can be produced in the final 



Running coupling corrections to inclusive gluon production 



3 



state would leave one factor of the coupling at an arbitrary renormalization scale 
jj,. Following [9] we rectify the problem by redefining the gluon production cross- 
section to include production of collinear gluon-gluon and quark-anti-quark pairs with 
the invariant mass lower than some collinear IR cutoff A 2 ^. This new observable is 
completely /^-independent and expressible in terms of the running coupling constants. 

3. Running coupling corrections to LO gluon production 

Gluon production at the lowest order in the coupling is shown in Fig. [1] in the A + = 
light-cone gauge. At this order the unintegrated gluon distribution is 

such that Eq. (pQ) reduces to 

da 2a 3 s C F 1 f d 2 q 

d?k T dy = 7T 2 fc 2 J q 2(k-q) 2 ' ^ ' 

Our goal is to set the scales for the three couplings in Eq. (jUJ). 




Figure 1. Diagrams contributing to the lowest-order gluon production in quark-quark 
scattering at high energy in the A + = light-cone gauge. 

The final result for the lowest-order gluon production cross section with the running 
coupling corrections included is [1] 

da _2C F a s (A 2 coll ) f d 2 q a 2 (q 2 ) a 2 ((fc - qf) 



d?k T dy 7T 2 k 2 J q 2 (k-q) 2 a s (Q 2 ) a s (Q* 2 ) ^ 
with the momentum scale Q being a complicated function of q and k — q given in pp. 
An interesting feature is that the scale Q 2 is complex-valued! The cross-section (I10p 
is, of course, real, as it contains a complex-valued coupling constant multiplied by its 
conjugate, a s (Q 2 ) a s (Q* 2 ). Eq. f JTU]) clearly looks like the fixed- coupling cross-section 
09]) with three factors of fixed- coupling replaced by the seven running couplings: we 
chose to refer to this structure as the septumvirate of couplings pp. 

An important feature of the scale Q is that in the q — > (or q — y k) limit 

i i fc2 1 

In -2- = In -2- + - (11) 
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reducing Eq. (ITO|) to 



da AC F a s (A 2 coll ) a s (k 2 ) a s (Q 2 s ) . k 2 



d 2 k T dy ~ 7T (tff Ql 



ln 7^- (12) 



4. Ansatz for the running coupling corrections in the &; T -factorization 
formula 

Using the lowest-order expression (fTOj) we would like to conjecture the following running- 
coupling generalization of Eq. ([I]): 



da 2C F 1 f 2 - - a s (A 



2 

coll J 



— / d*q <j>Jq, y) Mk -Q,Y-y) * v eM ' (13) 



d 2 k T dy 7T 2 fc 2 7 *^*""™v <*> yy a s (Q 2 ) a s (Q* 2 ) 
with the new (rescaled) distribution functions defined by (cf. Eqs. (J2J) and (jBJ) 



' A (k,y) = / d^re"*" V 2 iV G (r,b,y) (14) 

and 



= 7^3 | rf 2 W 2 re"^ V 2 n G (r,6 52/ ). (15) 

One may use this ansatz, along with Nq obtained from the rcBK evolution to further 
improve the existing CGC phenomenology for particle multiplicity along with its 
centrality and rapidity dependence at RHIC and LHC. 
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